Chiral gauge field and axial anomaly in a Weyl semi-metal 
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Weyl fermions are two-component chiral fermions in (3 + l)-dimensions. When coupled to a gauge 
field, the Weyl fermion is known to have an axial anomaly, which means the current conservation 
of the left-handed and right-handed Weyl fermions cannot be preserved separately. Recently, Weyl 
fermions have been proposed in condensed matter systems named as "Weyl semi-metals". In this 
paper we propose a Weyl semi-metal phase in magnetically doped topological insulators, and study 
the axial anomaly in this system. We propose that the magnetic fluctuation in this system plays the 
role of a "chiral gauge field" which minimally couples to the Weyl fermions with opposite charges for 
two chiralities. We study the anomaly equation of this sytem and discuss its physical consequences, 
including one-dimensional chiral modes in a ferromagnetic vortex line, and a novel plasmon-magnon 
coupling. 
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Introduction - In the quantum field theory, a (3+1)- 
dimensional massless Dirac fermion is decomposed to two 
independent two-component fermions known as the Weyl 
fermions. Weyl fermion has a definite chirality, left- 
handed or righthanded, determined by the sign of its 
spin polarization along the momentum direction. [l[ Clas- 
sically, the lefthanded and righthanded Weyl fermions 
are decoupled and can be coupled independently to two 
gauge fields, leading to a separate charge conservation. 
The gauge field that couples differently to Weyl fermions 
with two chiralities is called a chiral gauge field. For ex- 
ample the SU(2) gauge field in the Standard Model is a 
chiral gauge field. It is well-known that the chiral charge 
conservation is violated in a quantum theory of Weyl 
fermions in a background gauge field, which is known as 
the axial anomaly 0"3 ■ 

Recently, Weyl fermions are also introduced into con- 
densed matter physics. The Weyl fermions are shown 
to be the topologically robust boundary states of (4 + 
l)-d time- reversal (TR) invariant topological insulators 
(TI)[5j, and the axial anomaly corresponds to a topo- 
logical response of the (4 + l)-d TI. This approach is re- 
lated to the domain wall fermion approach^] and Callan- 
Harvey effect Q in high energy physics. By dimensional 
reduction, the (4+lV-d topological insulator is reduced to 
the (3 + l)-d TI@-|iq| and the Weyl fermion is reduced to 
(2 + l)-d surface states of the TI. Weyl fermions also ap- 
pear directly in (3 + l)-d gapless electron systems, which 
are named as "Weyl semi-metals" 0, Illl - i21j . Since a sys- 
tem with both TR and parity (P) symmetries have all en- 
ergy bands doubly degenerate, the Weyl semi-metal state 
can only be realized in a system breaking TR and/or P 
symmetry. 

A natural question is whether the chiral gauge field 
can be realized in the Weyl semi-metals, and if yes, what 
is the physical consequence. In this letter, we address 
these questions in TR breaking Weyl semi-metals. We 
show that generically a ferromagnetic moment couples 



to the Weyl fermions as a chiral gauge field. As an ex- 
plicit example system, we study a model of magnetically 
doped topological insulator, which can be driven into the 
Weyl semi-metal phase with strong enough magnetic mo- 
ments. The presence of the chiral gauge field leads to an 
anomaly equation satisfied by the charge current, which 
leads to new topological phenomena such as chiral one- 
dimensional states in a magnetic vortex, and a topologi- 
cal coupling between spin fluctuation and plasmons. 

Chiral gauge field and anomaly equation - We start 
with a general discussion of Weyl fermions in condensed 
matter physics. In a weakly interacting crystalline ma- 
terial, Weyl fermion states generically appear when two 
energy bands cross at a generic point Kq in the Brillouin 
zone. The low energy physics around Kq is described by a 
two-component Hamiltonian Hw — ^^2i j— x y z v ij^i cr ji 
with ki the momentum away from Kq, and Oj the Pauli 
matrices. The matrix V{j describes the generic linear 
coupling between momentum and spin degree of free- 
dom described by Uj. By rotating the basis one can 
always diagonalize Uy, and the three diagonal compo- 
nents are anisotropic velocities. Without losing general- 
ity, we restrict our discussion on isotropic Weyl fermions 
with the simple Hamiltonian H = hvfS-k. Our results 
on anomaly and chiral gauge field is insensitive to the 
anisotropy in the velocity. The sign of the Fermi veloc- 
ity Vf determines the chirality of Weyl fermion. A single 
Weyl fermion is topologically stable as long as translation 
symmetry is preserved (l2j. thus for a Weyl semi- metal 
with translation symmetry, any local perturbation can 
only move the nodal point Kq in the momentum space. 

According to the Nielsen- Ninomiya theorem (2^, [Hj|, 
in a lattice model the number of Weyl fermions with op- 
posite chiralities must be equal. Consequently, the min- 
imum number of Weyl fermions in a Brillioun zone is 2. 
Moreover, because TR symmetry preserves the chirality 
of Weyl fermion, in TR invariant system the minimum 
number of Weyl fermions is 4[l5j|. In the following, we 



focus on the "minimal Weyl semi-metal" which break TR 
but preserves P, with two Weyl fermions of opposite chi- 
ralities at wavevectors Kq and — Kq, related to each other 
by spatial inversion. Kq is a generic point in BZ away 
from TR invariant momenta. 

We consider an arbitrary perturbation to the system 
of two Weyl fermions. As long as the perturbation is 
so smooth that the momentum transfer is much smaller 
than 2|l£o|, the two Weyl fermions remains decoupled. 
The effective Hamiltonian of the lefthanded Weyl fermion 
under perturbation is Hl = tiVfa ■ k + SHl, with SHl 
a generic 2x2 Hermitian matrix. To the leading order 
one can ignore the k dependence and consider SHl as 
a constant term. Then SHl can always be expanded 
to the form SHl = hvfa ■ cll + a L with the last term 
proportional to identity. Adding this to the Weyl fermion 
Hamiltonian we find Hl = hvf{k + Sl) ■ a + ciql with 
a nL = (aoiiflL) behaving as a gauge field. Similarly 
one can define the gauge field a^R minimally coupled 
to the righthanded Weyl fermions with the Hamiltonian 
Hr = —hvf(k + an) ■ a + aoR. The two Weyl fermions 
can be described together by a 4 x 4 Hamiltonian: 



H = hvf [(k + A) ■ ar z + a ■ a) + a T z + A (1) 



with Af, = (&fj,L 



*Vr)/2 behaving like the electromag- 
netic gauge field, and a M = (a M £ — a fl n)/2 the chiral 
gauge field. and a M have different properties under P 
and TR. For example if the perturbation we consider is 
a fluctuation of a ferromagnetic moment, only a will be 
induced which is TR odd and P even. 

As known from the quantum field theory, when a Weyl 
fermion is coupled to a gauge field, the charge conser- 
vation is broken at the quantum field level, leading to 
the axial anomaly [H |24[, which can be described by the 



anomaly equation d^ L ^ = (-)^e^ f^ R) f, 
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Since the gauge field 
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where fpi = 9 M a !yL(fl) -d^^R) 
a nL(R) of lefthanded (righthanded) Weyl fermion is re- 
lated to the gauge field A^ and a M , the anomaly equation 



can also be rewritten as d fJ ,j' J ' L ^ R ' 
HhpXF^ + HUv), where 
tromagnetic field strength and f^ v 
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When A and a coexist, we find both the axial current 
and charge current are nonconserved with the anomaly 
equations 
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The equation @ is the axial current anomaly [3, [2l[ 
but with the additional term induced by chiral gauge 
field, while the equation Q indicates the conservation of 
charge current is also broken due to the combination ef- 
fect of chiral gauge field and electromagnetic field, which 
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FIG. 1: (a) "Chiral magnetic field" can be generated by 
the magnetic vortex configuration in a topological insulator 
cylinder. Here the vector b indicates the direction of the "chi- 
ral magentic field" . The Landau level spectrum of a massless 
Dirac fermion is plotted (b) for a uniform magnetic field B 
and (c) for a uniform "chiral magentic field" b. 



is the main focus of this paper. At the first glance, the 
breaking of the charge conservation seems something un- 
physical. Here we emphasize that the Weyl fermion de- 
scription is only a low energy effective theory and the 
high energy part is not taken into account. Let's de- 
fine j£ — —^ii >luXp a v d\A p , and the right-hand side of 
equation can be written as a total derivative of j£, 
and the charge conservation law <9 M (j M + j£) = is re- 
covered if 2b is regarded as a current from the high en- 
ergy part which is neglected in our description. Actu- 
ally we notice that the spatial component of jt~ is given 

by 3b = ^2^?® x E with the electric field E, exactly 
corresponding to the anomalous Hall response of Weyl 
fermion, as first derived in Ref 20]. To make our dis- 
cussion concrete, we first propose a realization of Weyl 
fermions and chiral gauge field in magnetically doped 
topological insulators, before discussing the physical con- 
sequence of this anomaly equation. 



Material realization - It is first suggested that Weyl 
fermions can be realized in pyrochlore iridatesfl2l|. and 
later another material HgCr2Se4 is also proposed [l3j. 
However both the materials include multiple Weyl 
fermions with the number larger than 2, making the sys- 
tem complicated, therefore it is desired to have a system 
with the minimal number of Weyl fermions, which actu- 
ally can be achieved by magentically doped topological 
insulators [lil [l6{. By substituting the atoms, it is pos- 
sible to tune the band gap of topological insulators, and 
even induce the phase transition between trivial and non- 
trivial phases, which has been realized in TlBi(Si_5Se5)2 
recently |25l - l27| . Near the transition point, the bulk gap 
is minimized and can be overcomed by the exchange cou- 
pling from magnetic doping. The ferromagnctism in the 
Cr or Fe doped Bi 2 Te3 and Sb2Te3 has been observed 
in experiment (281 - 130} . therefore the magnetically doped 
Bi 2 Se 3 and TlBiSe2 family of materials are the suitable 
platform for the realization of minimal number of Weyl 
fermions. Here we adopt the four band model[3ll |32| 
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FIG. 2: (a) The energy dispersion as a function of k z along 
a ferromagnetic vortex line in a topological insulators with 
the total angular momentum J z = |. (b) The corresponding 
radial wave function for the two chiral modes at k z = 0. Here 
red line is for the wave function near the vortex core (r=0) and 
the black line for the one at r=R. (c) The energy dispersion 
of the two zero modes for different J z = n + | . The gap is 
observed for large \n\, due to the finite size effect, which gives 
a cut-off of the total number of chiral modes. In (d), the wave 
function of the inner chiral mode moves outwards, hybridizing 
with the outer chiral mode (Here we take n = 3 and k z — 0.06 
1/A). The parameters of the four band model are taken to 



be M = 0, Mi = 0.342eV-A 2 , M 2 = 18.25eV- 



Bo 



1.33eV"-A, A = 2.82eV-A, U = O.leV and W = 0.06eV. 



with general mass terms, to describe these materials, 

H = H a + H 1 (4) 
Hq = e{k) + M(k)T 5 + B k z T 4 + A Q {k y Y 1 - k x T 2 ) 

ij 

where e k = C + Cik 2 + C 2 fcjj, M(k) = M + M x k 2 z + 
M2k 2 . The T matrices are defined as I\2,3 = o~x,y,zT X i 
T 4 = t v , T 5 = r z , and r Qb = [r a ,r h ]/2z (a, b = 
1,...,5). Ferromagnetism breaks T but preserves P, 
therefore by inspecting the symmetry property of T ma- 
trices (eg. the table III in the reference j32[), we im- 
mediately find only two sets of T matrices are allowed 
in Hi: = s ljk <Jk and T l4 = er^ (i,j,k = x,y,z). 
Generally T12 and can be induced by z-direction 
magnetization, while (r 14,^4) and (r23,r3i) originate 
from in-plane magnetization. It is shown that Ti4, T 2 4 
and Ti2 induce two Weyl fermions while T23, T3i and 
T34 yield a nodal ring [15]. Since now we are interested 
in the Weyl fermion regime, we focus on the simple 
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FIG. 3: (a) The experiment setup and the diagram for the 
corresponding physical process. Here the red line denotes 
the electromagnetic field, and the blue dashed line is for a z 
field, (b) The Landau levels for the four band model Q in 
the uniform magnetic field, (c) The imaginary part of the 
correlation function {a z (q)a z (—q)). Here we take u}q/u) v = 
0.1, v s = 0. 



case with Hi = UoT 12, yielding the energy dispersion 
E st = e(k) + s^jAKkl + kj) + (VM 2 + Blkl + t\U \) 2 , 
where s,t = ±1. The two bands in the middle with 
t = — 1 touches when the conditions M 2 + B^k 2 z = Uq 
and k x = k v = are satisfied. If we neglect the quadratic 
term in M for simplicity, the bulk gap is closed, realiz- 
ing Weyl f ermions, at the momentum k z — ±Kq with 
K = ^U 2 -M 2 if \U \ > \M \._ 

Next we consider the perturbation around the gap- 
less points k = (0, 0,±Ko) with the perturbed Hamil- 
tonian H' = BoSk z T y + Aq (6k y a x T x — 6k x a y T x ) + 
J2i= x ,y,z (/ i » cr * + v iVi T z): where p and v denote the 
magnetic fluctuation and 5k is the momentum ex- 
panded around (0,0, ±ifo)- We project the Hamil- 
tonian H' into the subspace expanded by |— ,— ) and 
|+,— ) by perturbation theory, obtaining the effec- 
tive Hamiltonian i? e // = hvf z {5k z a z T z + a z a z ) + 
h v f\\ J2i=x y {&ki<JiT z + a,i<Ti), where a denotes spin, r 
represents two Dirac cones at (0, 0, ±Kq), the Fermi ve- 
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Physical consequence - Let's introduce "chiral mag- 



netic field" b = V x a and "chiral electric field'' 



then equation ([3]) can be rewritten as 



e-B 



da 

dt ' 



(8) 



with p and j the charge density and current, respectively. 
In the following we study the physical consequences of the 
two terms on the righthand side of this equation. 
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The first term b ■ E describes the effect of a chiral 
magnetic field parallel to the electric field. The "chi- 
ral" magnetic field can also induce Landau levels, sim- 
ilar to the Landau levels by magnetic field. For a uni- 
form field b = b$e z , the Landau level spectrum is given 
by E± >a (ri) = zthvf^/k 2 + Ibtfn with n = 1,2,... and 
a = ± denote two Dirac cones. In addition, there 
are two zeroth Landau levels, both with the dispersion 
E a (0) = —hvfk z (a — ±), as shown in FigQ] (c). It 
should be noticed that the two zeroth Landau levels are 
one-dimensional modes with the same chirality, in con- 
trast to the case of an ordinary magnetic field B = B^e z 
shown in Fig. [T](b). In other words, the low energy dy- 
namics of the system is described by two chiral fermions 
in each area with one flux quanta of b. In this case, the 
anomaly equation dSj reduces to the chiral anomaly of 
ID chiral fermions |33L l34l| . 

A key difference of the chiral gauge field from the elec- 
tromagnetic gauge field is that the gauge vector potential 
a is physical and thus has to be single valued. A uniform 
b field corresponds to a a linearly increasing towards the 
boundary of the system, which is unphysical. Also the 
Landau level spectrum above does not consider the cou- 
pling between the two Weyl fermions. To obtain a more 
complete understanding to this problem, we consider the 
four band model Q with the magnetic vortex configura- 
tion mi4 = —Wo sin 9, 77124 = Wo cos 9 (9 is the angular 
coordinate), m\2 = Uq and all the other m<j = 0. As- 
suming a cylinder shape geometry as shown in Fig. (fTJ) 
(a), the system has translation symmetry in z direction 
and rotation symmetry according to z axis. Therefore 
the momentum k z and total angular momentum J z are 
good quantum numbers, and the Schordinger equation 
can be solved numerically for each sector of k z , J z by in- 
troducing a discretization, as described in the appendix 
in detail. Fig [5] (a) and (b) show the band dispersion 
along z direction for the angular momentum J z = 1/2, 
and the corresponding wave function. There are two gap- 
less modes with opposite chirality, different from the uni- 
form field Landau levels in Fig. [T] (c). The wavefunction 
shows that the two chiral modes are spatially separated, 
with one wave function around r — and the other one 
located at the boundary r = R, as shown in Fig [2] (a) 
and (b). From this result we see that the chiral modes in 
the zeroth Landau level are compensated by modes with 
opposite chirality on the boundary, which is expected 
since in such a finite system the number of left and right 
moving ID states must be equal. With increasing the an- 
gular momentum J z , the wavefunction of the inner chiral 
mode moves towards larger radius, as shown in Fig. [5] 
(d). For a finite system, the wave functions of the two 
chiral states near r — and r = R will overlap with each 
other for large angular quantum J z (Fig[2](d)), leading 
to a gap opening, as shown in Fig [2] (c) . Consequently 
we obtain a finite number of chiral modes in the zeroth 
Landau level. As expected, the number of chiral modes 
is determined by the total flux of a in the system, just 
like the case of uniform b field. In such a configuration, 



the consequence of the anomaly equation ([8j is actually 
a quantum Hall effect [20| . In an electric field E = Ez 
parallel to b, the anomaly equation describes a charge 
generation around the center of the system, while the 
charge on the boundary is annihilated. This is a conse- 
quence of a Hall current flowing along the radial direction 
towards the center, which can be measured in transport 
experiments. 

The second term on the right hand side of the anomaly 
equation © describes the combination effect of magnetic 
field and "chiral" electric field. To understand this term, 
consider a uniform magnetic field B = Bqz and a uni- 
form vector potential a = a z (t)z changing adiabatically 
in time. The anomaly equation leads to Sp — -^Sa z , 
with G = ^j^- the Landau level degeneracy. Therefore 
the change of a z leads to a charge density modulation 
proportional to it. To understand this equation we con- 
sider the Landau level spectrum for the four band model 
(|4|, as is shown in Fig[3](b). It is important to note that 
the two zeroth Landau levels with opposite chirality have 
the same spin polarization. Consequently, the exchange 
coupling of a = a z z with the zeroth Landau level states is 
equivalent to a scalar potential, which shifts the chemical 
potential and leads to the change of charge density. 

Since this term couples charge density and magnetiza- 
tion, it leads to an interesting physical consequence of the 
hybridization between the plasmon and magnon modes. 
The effective action of the present system can be given by 
Sa + S a + Srpa, where S A = J dkdajA (q)G A ^ j A (-q), 
S a = J dkduja z (q)G~ 1 a z (-q) and S RPA = 
J dkdoj [(A (q) + a z (q))U(A a (-q) + a z (-q))] with 
q = (co,k). Here a z has been rescaled to have the 
same dimension as Aq. Sa describes the dynamics 
of the scalar potential Aq with G^ 1 ~ v 2 k 2 with 
photon velocity v c , S a describes the dynamics of a z 
field with G" 1 ~ u 2 — v 2 k 2 — uj 2 where ujq gives the 
excitation gap and v s is the magnon velocity, and 
Srpa gives the effective action after integrating out the 
interacting fermions with the random phase approxi- 
mation, where n ~ — ^fc& 2 [35[ (u p is the plasmon 
frequency). Such a hybridization is shown schematically 
in Fig [3] (a). To see the effect on the dynamics of the 
magnon described by a Zl we integrate out Aq field and 

obtain S eff = J dkdua z (q) I G^ 1 + n+( ^-i ) a z (—q). 
The corresponding correlation function is given by 

uj 2 -v 2 k 2 ~uj 2 ~ Jll% +vq + iri) , 
which corresponds to the spin susceptibility. As plotted 
in Fig |3] (c), the correlation function has two poles, of 
which one corresponds to the intrinsic magnon excitation 
with the frequency around uq, while the other one only 
appears for finite k with the intensity proportional 
to k 2 and is induced by the plasmons with frequency 
around u) p . The plasmon frequency can be estimated as 
~ 35meV for Weyl fermions 36] with dielectric constant 
~ 100, Fermi velocity ~ 6.85 x 10 5 m/s, and electron 
density ~ 10 19 cm . Such an additional mode in 
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magnon spectrum can be observed in neutron scattering 
experiments and compared with the plasmon frequency 
Lo p determined by reflection spectroscopy. 
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Appendix A: Numerical method for the calculation 
of the chiral mode in the ferromagnetic vortex core 

In the appendix, we describe our numerical method for 
the calculation of energy dispersion and eigen wavefunc- 
tion for the ferromagnetic vortex configuration. We start 
from the four band model (4) in the main text and in the 
cylinder coordinate (r,6,z), the Hamiltonian takes the 
form of 



H = H + H 1 



(Al) 



H a = M{k)r z + B k z T y + A (k y a x T x - k x a y T x ) 



( M(k) 

M(k) 

iB k z iA k- 

V -iA k + iB k z 



—iBok z iAok- 
—iAok + ~iBok z 
-M(k) 
-M(k) J 



Hi = - Wo sin 9a x T z 

( U -iW e~ 
iW e ie -U 

\ 



Wo COS 9<7 y T z + Uq(7 z 
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-U J 
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-iW e lS 



-id. 



d, 
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where we have d x — cos9d r — ^r-dg and d y — s'm6d r + 
^de, therefore fc_ = k x - 
-i (cos9d r - s -^d e ) - (sm9d r 
^-d e , k + = -id x 
(sinddr + ^dg) = 



dr' 1 



1 a 



e -d e )- 
d„ = —i (cos 9d r 

+ c ^d e and k 2 x + k 2 y 
The above Hamiltonian has 



-ie id d, 



X 

-ie- iti d r 

•„2 , ;„2 



in-plane rotation symmetry along z axis and the cor- 
responding total angular momentum can be defined as 
J z = L z + \o z where L z = —i-§g and the Pauli matrix o z 
denotes the spin part. With the in-plane rotation sym- 
metry, the wavefunction ansatz can be taken as ip(r, 9) — 



[e'" 9 /i(r),e'("+ 1 ) 8 / 2 (r),e'" 8 / 3 (r), e 
the total angular momentum J. 
nian is changed to 



{n+1)e fi{r)] where 
n + | . The Hamilto- 



H = 



( M(n) + U -iW -iB k z 

iWo M{n + l)-U A (-d r + ^) 

iB k z A (d r + ^) -M(n) + U 

\A (-d r + f) iB k z -iW 



A (d r + ^) \ 
-iB k z 
iW 

-M(n + l)-U J 



(A2) 



where M(n) = M + Mik 2 z - M 2 + £) 

and the wave function is now given by ip = 
[fi ( r ) ! fz (r) , fz (r) , fi (r)] T . Let's introduce the new wave 



function ip as ip = -^tp, then the normalization relation 

/ rdrd9\ip\ 2 = 1 is changed to f drd9\ip\ 2 = 1, and the 
effective Hamiltonian is rewritten as 



H 



( M(n) + U -iW -iB k z Ao(d r + ^^-^\ 

iWo M(n+1)-U A (-d r + r -^) -iB k z 



iB k z A (d r + - 

\A (-d r + r -^) iB k z 



-1/2 



-M(n) + U 
-iW 



*W 

-M(n + 1)-U ) 



(A3) 



where M(n) = M + M\k 2 z - M 2 - . This 

Hamiltonian can be written in a compact form 
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H = 



Mq + Miki - M 2 



tP_ _ (n + 1/2) 1 



n + 1/2 

T z - M 2 5 d z T z 



+B Q k z T y + iA ^-a y T x + A ^— —a x T x + W a y T z + U Q a z 



We can discretize the Hamiltonian (|A3|) and solve 
the eigenstate problem for the radial equation numer- 
ically. The corresponding result is shown in Fig 2 of 
the main text. For n = 0, we indeed find two gap- 
less modes with the opposite velocities along z direc- 
tion, and these two gapless modes are spatially sepa- 
rated with one wave function mainly staying at r = 
and the other one at r — R, as shown by the red and 
black lines in Fig 2 (a) and (b) of the main text. How- 
ever with increasing n, a gap is opened between the two 
low energy modes, as shown in Fig 2 (c) in the main 
text. To get more analytical understanding of the radial 
equation, we consider the r — > oo limit with Uq = 0, 
k z = 0, where the radial Hamiltonian is simplified as 

H = (m - M 2 |^) t z + W a y T z + iAoa y T x £. With 

the wave function ansatz ip ~ e Xr cf>, we obtain the equa- 
tion A Q X4> = [(-Mo — M 2 A 2 ) a y T y + W t v ] 4> for the zero 
modes. Since \o y r y ,T y \ — 0, we can take the com- 
mon eigen-states of a y r y and r y for </>, a y T y 4>ts = t<f>ts 
and T y <j>t s = s4>t s , then the wave function can be ex- 
pressed as ip = Y, a ,t,s c a-ts e> "' i ' t,s)r< l ) t,s, with A given 
by X a (t,s) = -'^Wyt'W , The existence 



[1] M. E. Peskin and D. V. Schroeder, An Introduction 
to Quantum Field Theory (Addison- Wesley Publishing 
Company, Massachusetts, 1995). 

[2] S. L. Adler, Phys. Rev. 177, 2426 (1969). 

[3] J. Bell and R. Jackiw, Nuovo Cimento A 60, 47 (1969). 

[4] H. Nielsen and M. Ninomiya, Physics Letters B 130, 389 
(1983). 

[5] X.-L. Qi, T. L. Hughes, and S.-C. Zhang, Phys. Rev. B 

78, 195424 (2008). 
[6] D. B. Kaplan, Physics Letters B 288, 342 (1992). 
[7] C. G. Callan and J. A. Harvey, Nucl. Phys. B 250, 427 

(1985). 

[8] X.-L. Qi and S.-C. Zhang, Rev. Mod. Phys. 83, 1057 
(2011). 

[9] M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 
(2010). 

[10] J. E. Moore, Nature 464, 194 (2009). 

[11] G. Volovik, The Universe in a Helium Droplet (Claren- 
don Press, Oxford, 2003). 

[12] X. Wan, A. M. Turner, A. Vishwanath, and S. Y. 
Savrasov, Phys. Rev. B 83, 205101 (2011). 

[13] G. Xu, H. Weng, Z. Wang, X. Dai, and Z. Fang, Phys. 
Rev. Lett. 107, 186806 (2011). 

[14] A. A. Burkov and L. Balents, Phys. Rev. Lett. 107, 



(A4) 

I 

of the edge mode requires A+(+, +)A_(+, +) > or 
A + (+, — )A_(— , — ) > 0, leading to the following different 
regimes: in the normal regime MqM 2 > 0, the system has 
no zero mode when |Wo| < \Mq\ and one zero mode when 
| Ho I > I -Mo |, while in the inverted regime M M 2 < 0, the 
system has one zero mode when \Wo\ > \Mq\ and two 
zero modes when \ Wq\ < \Mq\. Taking into account the 
k z dependent term, it turns out that one zero mode case 
corresponds to the ID chiral state and two zero modes 
case is the ID helical state. However since time reversal 
is broken in the present system, the helical state is not 
protected and can be gapped. Therefore the only robust 
state is the chiral state when |Wo| > \M \. We emphasize 
that the transition at |Wo| = \Mq\ exactly corresponds 
to the condition for the appearance of the gapless Weyl 
fcrmions for the uniform magnetization. For the finite 
r, the terms proportional to £ and -j will push the chi- 
ral mode around r=0 outwards, thus with increasing the 
angular momentum number n, the wave function of the 
chiral mode near r = extends to the large r region and 
mixes with the chiral mode at r = R, opening a gap. 



127205 (2011). 

[15] A. A. Burkov, M. D. Hook, and L. Balents, Phys. Rev. 

B 84, 235126 (2011). 

[16] G. Y. Cho, e-print [arXiv:1110.1939l (2011k 

[17] C. Fang, M. J. Gilbert, X. Dai, and B. A. Bernevig, e- 

print larXiv:TTlT. 7309 (201 1). 
[18] J. Jiang, e-print larXiv:1112. 59431 (2011k 
[19] P. Hosur, S. A. Parameswaran, and A. Vishwanath, Phys. 

Rev. Lett. 108, 046602 (2012). 
[20] K.-Y. Yang, Y.-M. Lu, and Y. Ran, Phys. Rev. B 84, 

075129 (2011). 

[21] V. Aji, e-print laxXiv:1108.4426l (2011k 

[22] H. B. Nielsen and M. Ninomiya, Nucl. Phys. B 193, 173 

(1981). 

[23] H. B. Nielsen and M. Ninomiya, Nucl. Phys. B 185, 20 
(1981). 

[24] A. Zee, Quantum Field Theory in a Nutshell (Princeton 

University Press, New Jersey, 2010). 
[25] B. Yan, C. Liu, H. Zhang, C. Yam, X. Qi, T. Frauenheim, 

and S. Zhang, EPL (Europhysics Letters) 90, 37002 

(2010). 

[26] H. Lin, R. S. Markiewicz, L. A. Wray, L. Fu, M. Z. Hasan, 

and A. Bansil, Phys. Rev. Lett. 105, 036404 (2010). 
[27] S. Xu, Y. Xia, L. A. Wray, S. Jia, F. Meier, J. H. Dil, 



7 



J. Osterwalder, B. Slomski, A. Bansil, H. Lin, et al., Sci- 
ence 332, 560 (2011). 

[28] Y. L. Chen, J. H. Chu, J. G. Analytis, Z. K. Liu, 
K. Igarashi, H. H. Kuo, X. L. Qi, S. K. Mo, R. G. Moore, 
D. H. Lu, et al., Science 329, 659 (2010). 

[29] C. Chang, J. Zhang, M. Liu, Z. Zhang, X. Feng, K. Li, 
L. Wang, X. Chen, X. Dai, Z. Fang, et al., e-print 
larXiv: 1108.47541 (2011). 

[30] L. A. Wray, S. Xu, Y. Xia, D. Hsieh, A. V. Fedorov, Y. S. 
Hor, R. J. Cava, A. Bansil, H. Lin, and M. Z. Hasan, Nat 
Phys pp. 32-37 (2011). 

[31] H. Zhang, C. Liu, X. Qi, X. Dai, Z. Fang, and S. Zhang, 
Nat Phys 5, 438 (2009). 



[32] C.-X. Liu, X.-L. Qi, H. Zhang, X. Dai, Z. Fang, and S.-C. 
Zhang, Phys. Rev. B 82, 045122 (2010). 

[33] K. Johnson, Physics Letters 5, 253 (1963). 

[34] R. Jackiw, Topological investigations of quantized gauge 
theories, in Current Algebra and Anomalies, edited by 
Sam B. Treiman, Roman Jackiw, Bruno Zumino and Ed- 
ward Witten, Princeton University Press, 1985. 

[35] N. Nagaosa, Quantum Field Theory in Condensed Matter 
Phystcs (Springer, Berlin, 1999). 

[36] S. Das Sarma and E. H. Hwang, Phys. Rev. Lett. 102, 
206412 (2009). 



